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Abstract 



A Monte Carlo algorithm for computing quantum mechanical expectation 
values of coordinate operators in many body problems is presented. The algo- 
rithm, that relies on the forward walking method, fits naturally in a Green's 
Function Monte Carlo calculation, i.e., it does not require side walks or a bilin- 
ear sampling method. Our method evidences stability regions large enough to 
accurately sample unbiased pure expectation values. The proposed algorithm 
yields accurate results when it is applied to test problems as the hydrogen 
atom and the hydrogen molecule. An excellent description of several prop- 
erties of a fully many body problem as liquid ^He at zero temperature is 
achieved. 
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I. INTRODUCTION 



Quantum Monte Carlo (QMC) methods have become an invaluable tool in the study 
of many-body systems over the last decades. Among them, the Green's Function Monte 
Carlo (GFMC) methodEHI has been extensively applied to the calculation of ground-state 
properties of small molecules and quantum liquids and solids at zero temperature. Within 
the GFMC techniques one can distinguish between the domain GFMC,0 which stochastically 
constructs the Green's function, and the Diffusion Monte Carlo (DMC) method^ based on 
a short-time approximation for the Green's function. We will focus our discussion on the 
DMC method but the algorithm we present here for the evaluation of pure estimators can 
also be easily incorporated in a domain GFMC program. 

The DMC method solves the Schrodinger equation in imaginary time for the function 
/(R,t) = ^(R)vl/(R,t), 

- = -\ V^/(R, t) + \ Vr(F(R)/(R, t)) + (E,(R) - E) /(R, t) , (1) 

being \1/(R, t) the wave function of the system and a trial function used for importance 
sampling. In Eq. (|l]), which is written in atomic units. El = ipCR.)'^ HipCR) is the local 
energy and F(R) = 2'?/'(R)~-'^Vr-?/'(R) is the so-called quantum force; R stands for a 3A^- 
coordinate vector and E is an arbitrary energy shift. The Schrodinger equation for /(R, t) 
(HI) presents in the right-hand side three terms that are associated, by analogy to classical 
equations, to diffusion, drift and branching processes, respectively. The asymptotic solution 
of Eq. ([H), for any value E close to the energy of the ground state and for long times 
{t —>■ oo), gives the ground-state wave function $0(1^) provided that there is a nonzero 
overlap between \1/(R, t = 0) and $o(R)- The formal solution of Eq. (|l]) is 

/(R', t + At) = JdR G{R', R, Ai)/(R, t) , (2) 

where the Green's function G(R', R, At) gives the probability of transition from R to R' 
in a time interval At. The DMC method solves stochastically Eq. (Q) assuming reasonable 



2 



approximations for the Green's function when At O.uO After an iterative process, the 
asymptotic solution /(R, t oo) = '?/'(R)$o(R') is finally obtained. 

The direct calculation of the expectation value of an operator v4(R) from the asymptotic 
function /(R, t —>■ oo) corresponds to the mixed estimator 

WR)|A(R)|$o(R)) 
= (^(R) I $o(R)) • 

It gives an exact result only when A is the Hamiltonian H or commutes with H. Among 
the different methods to calculate expectation values for operators that do not commute 
with H, the extrapolation method^'! is the most widely used. Following this method, which 
has been extensively applied in QMC calculations, one has an approximation to the "pure" 
(exact) value, 

($o(R)|A(R)|$o(R)) 

^^^^^^^ - ($o(R) I $o(R)) ' 

by means of a linear extrapolation 

(A(R))e = 2 (A(R))^ - (A(R)), , (5) 

where 

(^(R)|A(R)|V;(R)) 

(^(R))^ = / ,/-DM , (6) 



is the variational estimator of yl(R). 

The accuracy of the extrapolation method is closely related to the trial wave function used 
for importance sampling. Furthermore, in spite of using accurate trial wave functions, the 
extrapolated estimator is always biased in a quantity difficult to assess. In order to overcome 
these important restrictions, several algorithms have been proposed in the last years. In the 
approach of Zhang and Kalofi a bilinear sampling method is used. In this scheme, the 
system is doubled and the random walks take place in an enlarged configuration space. 
Other approaches are based on the estimation of the quotient ($o/ 4') from the asymptotic 
offspring coming from the branching term.i In this line, Barnett et a/.00 and Runge and 



Rungetil have constructed tagging algorithms to properly account for the asymptotic number 
of descendants. Although the later scheme has provided satisfactory results in some specific 
cases, the large fluctuations observed in the asymptotic offsprings, and therefore in the 
corresponding weights ($0/ '^), have precluded the consideration of the forward walking as 
a stable and reliable method. In contrast with these considerations, we find that these 
statistical fluctuations (of unphysical origin) show a highly depressed effect over integrated 
quantities, and that in order to accurately sample pure estimators stable regions can be 
reached. The method we present is somehow related to the one of Ref. |] but with the 
advantage of not requiring a tagging algorithm. As we shall show in Section 2, the averages 
are basically taken as mixed expectation values, and therefore the pure estimators can be 
readily incorporated in the original Monte Carlo algorithms. 

The layout of this paper is as follows. In Section 2 the algorithm for the evaluation of 
pure estimators of coordinate operators is described. In order to check the correctness of 
our implementation in a Diffusion Monte Carlo code, as well as the capability and resolution 
of the method, results for several moments of H and II2, where exact results are available, 
are presented in Section 3. Pure results for the partial energies and structure properties of 
liquid ^He are reported in Section 4. The application of the method to a real many body 
problem is a compelling test for its reliability and stability. Finally, the main conclusions 
are discussed in Section 5. 



II. PURE EXPECTATION VALUES 



The pure estimator of an operator A(R) (|) may be written as 



(A(R))p= <|.o(R) 



A{R) 



^(R) 



^(R) 



<fo(R) 



'^o(R) 
V^(R) 



V^(R) 



(7) 



Following Liu et a/.,i $o(R)/^(R) can be obtained from the asymptotic offspring of the R 
walker. In fact, assigning to each walker Rj a weight ly(Rj) proportional to its number of 
descendants 
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W(R) = n(R,t^ oo) , (8) 

Eq. (0) turns out to 

where the summatory J2i runs over all walkers and all times in the asymptotic regime. As 
it is clear from its proper definition, the weight of a walker existing at time t, given by Eq. 
(1), is not known until a future time t' > t + T, being T a time interval long enough so that 
Eq. (D could be replaced by W{K{t)) = n{R{t')). 

In order to proceed to the evaluation of Eq. @ two different approaches are possible. 
In the first one, a tagging algorithm capable of identifying, at any time, which walker of 
any precedent configuration originated an actual walker could be used. Then, one could 
determine the number of descendants of the former Rj, and accumulate its contribution to 
Eq. (|^) "from the distance" . Such a tagging algorithm has been devised in Refs. |9|-pn]. On 



the other hand, one can work out an algorithm that operates with only the actual values 
of A(Ri), in such a way that a weight proportional to its future progeny is automatically 
introduced. This second approach is the one we have followed in the present work. 

The schedule of the algorithm is the following. The set of walkers at a given time {Rj} 
and the values that the operator A takes on them {Ai} evolve, after a time step, to 

m - m (10) 
{A} - {a:} . (11) 

In the same time interval, the number of walkers changes to A^'. In order to sample the 
pure estimator of A, we introduce an auxiliary variable {Pi}, associated to each walker, with 
an evolution law given by 

{P^} - {PI} = m + {P-} , (12) 
where {-P/} is the old set {Pi} "transported" to the new one, in the sense that each element 



Pi is replicated as many times as the Rj walker, without any other changes. {Pi} is initialized 
to zero when the run starts. 



With this procedure, after M addition steps (p^ ) we end up with a set of Nf values {Pi}- 
A pure estimator of A is given by 

Nf 

{A{K)), = Y.{P.}KMxNf) . (13) 

i=l 

The contribution to the {Pi} entering in Eq. (|l^), coming from a generic walker R(t), 
can be determined following the evolution of the series. Clearly there is no contribution 
from R(t) until time t is reached. At this moment, A(R(t)) enters in one of the rows of 
{P} (p!2D. From now on, if any of the descendants of R(t) disappears or replicates, the 
former contribution so does. As a result, y4(R(t)) appears in as many rows of {P} as 
descendants of R(t) exist, and therefore its contribution to Eq. (|^) is proportional to the 
weight W(R{t)){^. 

A final regard concerning the implementation of the algorithm has to be made. In Eq. 



(|T^), the "transport" operation accounts for the replication of the v4(R) contribution. In 
order to ensure the asymptotic condition (^, the series are continued for a while only with 
the reweighting law 

{p^} - {PI} = {pn ■ (14) 

Since a calculation is usually divided in several blocks, one can collect data during a block and 
allow for a further reweighting in the following one. In this second block, new information 
can be accumulated to be reweighted in the next block. This mechanism can be incorporated 
in the algorithm in a rather simple way. The final result is that, after a first initialization 
block, each new block gives a value for the pure expectation value of A. 

An alternative to the simple branching algorithm, implicitly assumed in the above 
method, is the use of weights p(Rj) related to the branching factor. In fact, it has been 
provedi'0 that the branching algorithms with weighting allow for some reduction in the 
variance of the expectation values. Our method for computing pure estimators is easily 
extended to these algorithms. In particular, the evolution laws become: 
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{P^} - {PI} = X A[} + 



X P* 



{^.} - {PI} 



X P* 



(15) 
(16) 



whereas the expression of the pure expectation value (|T^) is only modified by a normalization 
factor. 



III. APPLICATION TO SIMPLE SYSTEMS: H AND H2 

As a test of the algorithm developed in the preceding Section, we present results for the 
lowest coordinate moments of the hydrogen atom H and the hydrogen molecule H2. In both 
systems the nuclei are kept fixed and relativistic corrections are neglected. Atomic units 
have been used in all the Section. 

The DMC program used in the calculations is exact up to order (At)^ in the short-time 
approximation for the Green's function. More specific details of the algorithm are given 
m Ref. O Chini has extensively discussed the relation between a quadratic time-step 



dependence of the eigenvalue and the violation of the cusp condition in electronic systems. 
However, the achievement of a quadratic dependence in At, which has also been discussed 
by Umrigar et a/., Ill is not the main objective of the present work. Our aim is to check the 
efficiency of the algorithm for the extraction of pure expectation values in simple systems, 
as H and H2, where exact results are available. 

We have tested the reliability of our method including it in two versions of the original 
code, corresponding to the use or not of the weights p(Ri) (15,16). In both cases, satisfactory 
results are obtained, the variance of the mixed and pure estimators being slightly reduced 
when the weights p(R.i) are considered in the branching process. 
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A. Hydrogen atom 



Two different guiding functions are used for importance sampling in the H calculation. 
The first one corresponds to a Is Slater orbital 



with a slightly modified exponent {a = 0.9). The second one is taken as the product of the 
Is Slater orbital by a gaussian 



with a = 1.0 and (3 = 0.06. The analytic variational energies are Ej; = —0.495 for ipj 
and = —0.4853 for ipjj, to be compared with the exact result E = —0.5. A difference 
between ipj and ipii, which could be relevant in the time-step dependence of the energy, is 
that whereas ipu satisfies the cusp condition ipj does not. As far as the DMC calculation 
is concerned, we have used At = 0.05 in both cases with no significant differences when 
At is reduced by a factor of two. The number of walkers was maintained at a value of 700 
with an unnoticeable bias respect to larger populations. The samplings were performed over 
approximately 4 ■ 10^ configurations. 

In Table I variational (^), mixed @ and extrapolated estimators of the potential 
energy V, the radial distance r, the squared radial distance and z"^ are reported in com- 
parison with the exact results. The extrapolated expectation values improve the mixed 
results lying near the exact ones. However, some differences which depend on the trial func- 
tion used for importance sampling remain, showing that the extrapolation method suffers 
from a systematic bias related to ip. 

The pure expectation values of the same coordinate moments are reported in Table II. 
Neither the ipj nor the ipu results are biased with respect to the exact values. In fact, as it 
happens in the exact mixed estimator for the Hamiltonian, the quality of the trial function 
is only reflected in the magnitude of the variance. This influence may be observed in the 
larger errors of the pure estimators for tpjj with respect to the ones for ipj. 



tpi = exp(— ar) 



(17) 



■0// = exp(— ar — /3r^) 



(18) 
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The DMC calculation is divided in blocks of a number of iterations AL. According to 
the algorithm developed in Section 2 the block length has to be long enough to ensure the 
pure estimation in the asymptotic regime. In Figs. 1 and 2, the AL-dependence of the pure 
expectation value for is plotted for ipi and ipu-, respectively. Also shown are the exact 
result (solid line) and the extrapolated estimator (dashed line) corresponding to the trial 
function used in the calculation. The bias coming from the wave function components other 
than the ground state is rapidly suppressed, as expected from the evolution law exp(— i/t). 
The asymptotic condition is satisfied in both cases for values AL > 500. Beyond a transient 
regime, the prediction of the pure estimator is stable for a wide range of AL values with 
a negligible systematic bias. The statistical error in the ipn case is larger than in the ipi 
one but, in both cases, the central value reproduces accurately the exact results. On the 
other hand, the extrapolated predictions are biased respect to the exact and pure values, 
significantly for ipu as expected from its poorer variational quality. 



B. Hydrogen molecule 

The trial wave function we have used in the study of the hydrogen molecule is of the 
formS 

^ = (piriA, Tib) 0(r2A, ^25) exp(ari2/ (1 + feria)) (19) 

with the molecular orbital 

0(nA, TiB) = exp(-C riA) + exp(-C n^) . (20) 

The distances riA, ^iB correspond to the electron-nucleus separation, and stands for 
the electron-electron distance. The internuclear separation is kept fixed in the equilibrium 
distance r^B = 1-401. In order to test the accuracy of the pure algorithm, we have used the 
trial function (p!9|-pO|) with two different sets of parameters:@ 

= $(( = 1.189, a = 0.50, b = 0.40) (21) 



and 



= $(( = 1.189, a = 0.0, b = 0.0) . (22) 

The value of ( is obtained from the cusp condition between an electron and a nucleus {( = 
1 + exp(— r^^)). The Jastrow factor between the electrons appearing in the general form 
(^) is suppressed in ipu whereas it is considered in ipj with a value for the parameter a which 
guarantees the electronic cusp condition. The variational energies are E^^ = —1.1471(9) for 
and = —1.1288(8) for ipn, to be compared with the exact result E = —1.17447. . ..0 
The DMC calculations have been carried out with At = 0.05 and an average population 
of 500 walkers. The sampling has been made over 2 ■ 10^ configurations. We have checked 
that the biases due to the time step and the finite population are smaller than the statistical 
error. 

The H2 results closely follow the trends observed in the H calculation. The AL- 
dependence in the pure estimator of is shown in Fig. 3 for ipj, and in Fig. 4 for ipjj. 
Similar behaviours have been obtained for the other operators. One can see that the asymp- 
totic region is already reached at relatively small values of AL. In the asymptotic regime, 
the pure values fluctuate around the exact valued (solid line) with a statistical noise related 
to the quality of the trial function used for importance sampling. The growth of the error 
bars due to the death of walkers in the forward walking process is only significant for the 
largest AL values. The extrapolated prediction (dashed line) is manifestly biased. 

In Table III, we report results for variational, mixed and extrapolated estimators of V, 
and z"^ using ipi ( pT]) and ipu (^2]) as trial functions. As one can see, the variational results 



of the coordinate moments for ipu are closer to the exact values than for ipj, although ipi 
is energetically preferred to ipji. Nevertheless, the simple extrapolated expectation values 
are in neither case statistically compatible with the exact values.^ In Table IV the pure 
expectation values for the same moments are reported in comparison with the exact results. 
In spite of the shortcomings of these trial functions in reproducing the properties of H2, the 
pure estimator does reproduce the exact values independently of the trial wave function. 
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IV. LIQUID ^HE 



Domain GFMC and DMC have been extensively applied to the study of the ground state 
properties of liquid '^He for the last fifteen years.00 The exactness of these methods joined 
with the accuracy in the knowledge of the He interatomic potential have made feasible an 
excellent agreement between theoretical results and experimental data. In order to sample 
expectation values other than the Hamiltonian, e.g., partial energies or the radial distribution 
function g{r), the extrapolation method has been commonly used. In spite of the success in 
describing properties as g{r) a small bias is present in integrated quantities as the potential 
energy. Furthermore, the extrapolated estimator has evidenced its shortcomings in the 
calculation of density profiles of ^He clusters producing unphysical negative values for p(r) 
in the cluster surface.0 We have applied the algorithm developed in Section 2 to bulk liquid 
^He in order to show both its applicability to a fully many body problem and its capacity 
of removing the uncertainties introduced by the extrapolation method. 

The Schrodinger equation is solved by means of a Quadratic Diffusion Monte Carlo 
method considering the A^-particle Hamiltonian 

H = -^Vl + V{R) , (23) 

being R = (ri, . . . ,r]\f) and V(R) the interatomic potential. The results presented below 
have been obtained considering the HFD-B(HE) potential proposed by Aziz et a/.0 In a pre- 
vious paper0'Ei we have shown that this renewed version of the well-known Aziz potential^ 
improves appreciably the equation of state of liquid "^He with respect to the Aziz results. 

In order to establish the influence of the trial wave function used for importance sam- 
pling several options have been considered. The first one is the simple McMillan two-body 
function^ 

51 



ipji = n 



11 

2 Ir 



(24) 



with b = 1.20 0" (cr = 2.556 A). The second one is an improved version of ( ^^ proposed by 
Reatto@ 
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ipj2 = n 

i<j 



2 \r 



L 



exp 



r^3 - A 



(25) 



with b = 1.20 a, L = 0.2, A = 2.0 a, and A = 0.6 a. Finally, we have also used a trial wave 
function which contains three-body correlationsil 



ipjT = V'Ji exp 



where 



k i<j 



l^k 



(26) 



(27) 



and 



^(r) = exp 



r — Tf 



(28) 



The values for the triplet parameters are A = —1.08 cr~^, = 0.80 a and = 0.41 a. In the 
three trial functions, the values of the parameters are those which optimize the variational 
energy at the experimental equilibrium density Po^^' = 0.365 a^^. All the results pre- 
sented below correspond to the density p^^^ that coincides with the theoretical equilibrium 
density.El 

Results for the potential energy per particle using ipji, ipj2 and ipJT as importance 
sampling are reported in Table V. A small but significant difference between the extrapolated 
results {V/N)e appear, pointing to a bias related to the quality of the trial wave function. 
The bias is completely removed when the pure estimator is calculated, as one can see in the 
last row of Table V. The three values for {V/N)p are indistinguishable and, what is more 
important, they evidence a systematic error of the extrapolation approximation. In fact, 
none of the extrapolated values is statistically compatible with the common pure value, being 
the closest estimation the one obtained with ipjT which actually is the best variational choice. 
Considering the result for the energy per particle E/N = — 7.267±0.013 K,liithe pure result 
for the kinetic energy is T/N = 14.32 ± 0.05 K. Experimental determinations from analysis 
of deep inelastic scattering data predicts a slightly lower value (T/A^)'^^^* = 13.3 ± 1.3 
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being the difference mainly due to the significant errors in the experimental measurement 
of the tail of the response function. 

As far as the stability of the method is concerned, the dependence of the pure expecta- 
tion value of V/N on the length of the forward walking is plotted in Fig. 5. The results 
obtained (points with errorbars) follow the trends observed in Section 3 (Figs. 1-4). After 
a transition regime, and already for relatively small AL values (AL > 250), an asymptotic 
limit is reached where the systematic error is practically negligible. Notice that in the simple 
algorithm we have presented in Section 2 a forward walking of length AL is constructed from 
data ranging from L to 2L, and hence the length of the forward walking is not the same 
for all the walkers. This effect is not relevant provided that a region of stability exists. On 
the other hand, one can determine the asymptotic value within a single run collecting data 
for different block lengths. The statistical errors associated to each individual AL value 
can be lowered in the same way as mixed estimators, i.e., continuing the evolution of the 
series. The biases associated to the extrapolated expectation values are also shown in Fig. 
5, where {V/N)^ using '(pj2 and ipjx are represented by a long-dashed and a short-dashed 
lines, respectively. 

Other important quantities in the study of quantum liquids can also be calculated with 
the pure algorithm. In particular, the two-body radial distribution function 



NiN-1) /|<l>o(ri,...,rjv)Nr3...rfr 



/ |$o(ri, . . . ,rAr)|^rfri . . .rfrTv 

and the static structure function 



with 

N 

Pq = Ee^''-^^- (31) 

i=l 

The result obtained for g{r) is shown in Fig. 6 in comparison with the experimental data of 



Ref. |2^. As one can see, the pure expectation value of g{r) is in a good agreement with the 
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experimental g{r) for all the calculated r values. In Fig. 7 the pure structure function S{q) 
is plotted together with the experimental measures of Refs. p5| , |26| . An overall agreement 
between the theoretical and experimental S{q) is obtained, lying our result well between 
the two experimental determinations. The extrapolated estimations of g{r) and S'(g)Eilii 
are not significantly different of the pure result. It is clear that the difference between the 
results provided by the extrapolated and pure estimators is larger for integrated quantities 
as, for instance, the partial energies. 



V. CONCLUDING REMARKS 

We have presented a simplified version of the forward walking method to obtain unbiased 
expectation values for operators that do not commute with the Hamiltonian. One of the 
advantages of this algorithm, in front of others based on side walks or bilinear sampling, 
is that it enters "naturally" in any domain or short-time Green's function Monte Carlo 
program. The sampling of the pure estimator closely follows the standard procedure to 
sample mixed expectation values. Nevertheless, the main point is the stability of the method. 
The exponential decrease of the bias with the forward walking length and the evidence of a 
slow increase of the statistical uncertainties for physical quantities result in stability regions 
large enough to sample pure expectation values with negligible biases. In contrast, we do 
find also the common result that the weight associated with the offspring of an individual 
walker fluctuates all along the stability region. However, looking for stability in the values of 
the weights is more than what should be asked for. What is mainly required from a Monte 
Carlo method is stability in the expectation value of an operator A, {A)p. The computation 
of this average has clearly a much better chance of success. 

The accuracy of the method has been first verified in the H atom and the H2 molecule. In 
both systems, the pure expectation values reproduce the exact results with statistical errors 
which are not appreciably larger than the ones associated to the extrapolated predictions. 
In all cases, the extrapolated expectation values appear significantly biased respect to the 
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pure/exact values in a quantity which is related to the trial function but difficult to assert 
a priori. 

Finally, the pure algorithm has been applied to study some properties of liquid ''He at 
zero temperature. The implementation of this algorithm in a many body problem is also 
quite straightforward and the results obtained follow the same trends analysed in the simple 
systems (H,H2). The method is stable and generates results which are not biased by the 
importance sampling as it happens with the extrapolated estimations. In order to reduce 
the error bars of the pure values to the level of those associated to the total energy, the series 
have to be a bit longer. However, one does not have to perform the auxiliary Variational 
Monte Carlo calculation required by the extrapolation method and, more importantly, the 
guarantee of an exact result is fulfilled. 
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FIGURES 

FIG. 1. Pure expectation value of for H as a function of the block length AL using ^/^/. The 
solid and dashed lines correspond to the exact and extrapolated results, respectively. 



FIG. 2. Same as in Fig. 2 but for ipjj. 

FIG. 3. Pure expectation value of for H2 as a function of the block length AL using ipj. The 
solid and dashed lines correspond to the exact and extrapolated results, respectively. 



FIG. 4. Same as in Fig. 4 but for ^/j. 

FIG. 5. Pure expectation value of V/N for liquid ^He at p'^^^ as a function of the block length 
AL. The long-dashed and short-dashed lines stand for the extrapolated estimations using ijjj2 
and ipjT, respectively. 

FIG. 6. Pure expectation result of the two-body radial distribution function (solid line) for 
liquid ^He at Pq^^^ in comparison with the experimental data (points) of Ref. 25. 



FIG. 7. Pure estimation of the static structure function (points) for liquid ^He at /Oq^^* in 
comparison with the experimental determinations of Refs. ^ (solid line) and ^ (dashed line). The 
error bars of the theoretical points are only depicted where larger than the size of the symbols. 
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TABLES 

TABLE I. Variational, mixed and extrapolated expectation values for H using ■^j and ■0//. All 
the results are analytic. 
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TABLE II. Pure expectation values for H, using ^/ and V'//) in comparison with the exact 
results. The statistical errors are indicated in parentheses. 

00 (r) (r^) (^^) 

i)! -0.9987(10) 1.4999(10) 3.0002(36) 1.0004(17) 

-0.9975(14) 1.4993(28) 2.995(14) 0.9964(61) 

Exact -1.0 1.5 3.0 1.0 
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TABLE III. Variational, mixed and extrapolated expectation values for H2 using ipi and ipu. 



{V) 



Variational 
Mixed 

Extrapolated 
tpn 

Variational 
Mixed 

Extrapolated 
Exact 



-2.1034(9) 
-2.2415(7) 
-2.3796(17) 

-2.2254(19) 
-2.3012(13) 
-2.3770(32) 
-2.3489 



3.0740(1) 

2.7530(14) 
2.4320(28) 



1.2813(38) 
1.1281(12) 
0.9749(45) 



2.6155(23) 
2.5416(42) 
2.4677(87) 
2.5464 



1.0942(22) 
1.0409(19) 
0.9876(44) 
1.0230 



TABLE IV. Pure expectation values for H2, using tpi and ipn, in comparison with the exact 

results. 

{V) {r^ jz') 

tpi -2.3448(24) 2.5424(46) 1.0244(23) 

ipii -2.3454(39) 2.5412(74) 1.0210(44) 

Exact -2.3489 2.5464 1.0230 



TABLE V. Variational, mixed, extrapolated and pure expectation values of V/N (in K) for 
liquid ^He at p^"^^ using different trial wave functions ip. 

IpJl "072 IpJT 

{V/N)y -21.054(26) -21.311(18) -21.348(20) 

{V/N)m -21.459(8) -21.600(8) -21.541(8) 

{V/N)e -21.864(30) -21.889(24) -21.734(25) 

{V/N)p -21.56(5) -21.59(5) -21.58(5) 
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